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Abstract 

On anisotropic lattices with the anisotropy ^ = aa/ttr the follow- 
ing basic parameters are calculated by perturbative method: (1) the 
renormalization of the gauge coupling in spatial and temporal direc- 
tions, ga and gr, (2) the A parameter, (3) the ratio of the renormalized 
and bare anisotropy 77 — ^/^b and (4) the derivatives of the coupling 
constants with respect to ^, dg~^/dS, and dg~^/dS,. We employ the 
improved gauge actions which consist of plaquette and six-link rectan- 
gular loops, co-P(l X 1)^^ -I- ciP(l X 2)^j,. This class of actions covers 
Symanzik, Iwasaki and DBW2 actions. The ratio 77 shows an impres- 
sive behavior as a function of ci, i.e., 77 > 1 for the standard Wilson and 
Symanzik actions, while 77 < 1 for Iwasaki and DBW2 actions. This 
is confirmed non-perturbatively by numerical simulations in weak cou- 
pling regions. The derivatives dg~^/d£, and dg~^/d£, also change sign 
as — ci increases. For Iwasaki and DBW2 actions they become oppo- 
site sign to those for standard and Symanzik actions. However, their 
sum is independent of the type of actions due to Karsch's sum rule. 
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1 Introduction 



Anisotropic lattices allow us to carry out numerical simulations with the fine 
temporal resolution while keeping the spatial lattice spacing coarse, i.e. Ur < 
Cct, where aj- and are lattice spacing in the temporal and spatial direc- 
tions, respectively. This is especially important for QCD Monte Carlo sim- 
ulation at finite temperature and heavy particle spectroscopy. There have 
been many such calculations like glue thermodynamics Q, hadron masses 
at the finite temperature[^, glueballsHI and heavy quark spectra[^. The 
anisotropic lattice may become an important tool for the calculation of the 
transport coefficients of the quark gluon plasmafs], § and for the deter- 
mination of spectral functions at finite temperature 0, |8|. In numerical 
simulations on the anisotropic lattice we need the information upon the 
renormalization of anisotropy which is given by r] = £,/£,b where ^ is the 
renormalized anisotropy, a^/ar, and is the bare one. Karsch has first 
studied anisotropic lattices perturbatively with the standard plaquette ac- 
tion and obtained r] together with anisotropy coefficients which are defined 
by the derivatives of the spatial and temporal gauge couplings with respect 
to i.e., dg~'^/d^ and dg~'^/d^, and the QCD A parameter p. In Ref.[|lO| r] 
was determined non-perturbatively by analyzing Wilson loops in numerical 
simulations. 

Anisotropic lattices play an essential role in the analysis of thermody- 
namics of QCD. To get the thermodynamics quantities like the internal 
energy and the pressure from numerical simulations, one needs to know the 
anisotropy coefficients. 

On the anisotropic lattice, we can change the temperature by changing 
not only A'',- but also a,-. This allows us to adjust the temperature continu- 
ously with fixed spatial volume. 

Recently it has been recognized that improved actions are effective to 
get reliable results from lattice QCD simulations on the relatively coarse lat- 
tice; lattice artifacts due to the discretization are expected to be much less. 
Therefore, it is important to employ improved actions in the anisotropic lat- 
tice calculations, since the lattice is rather coarse in spatial direction in cur- 
rent numerical simulations. Garcia Perez and van Baal pursued first this di- 
rection, i.e., they have determined the one-loop correction to the anisotropy 
for the square Symanzik improved action [ 11 1. 

In this article we study the improved actions which consist of plaquette 
and six-link rectangular loops, 

5 (X ^ [coP(l X 1)^, + ciP(l X 2)^,] , (1) 
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where cq and ci satisfy the relation cq + 8ci = 1. This class of actions 
covers tree level Symanzik action without the tadpole improvement (ci = 
|l^ Iwasaki action (ci = -0.331) Q and DBW2 action (ci = 
—1.4088) |jl5|. They are most widely used in simulations of recent days. For 
the class of actions which consists of planar loops the anisotropic lattice can 
be formulated in the same way as for the standard plaquette action. This 
may not be the case for improved actions which include non-planar loops in 
three or four dimensions. 

In the following, we will calculate the A parameter, 77 = ^/^b and the 
anisotropy coefficients, in weak coupling regions mainly by perturbative cal- 
culations. 

In sect. 2 we briefly review the formulation of the anisotropic lattice 
with the improved actions and summarize formulae which will be used in 
this paper. In sect. 3, we outline the background field method and discuss 
the removal of the infrared divergence. In sect. 4 we present results of the 
perturbative calculations. The ci dependence of r/ and anisotropy coefficients 
are studied in detail. Since the behavior of rj is very important for practical 
use, we will study it by the numerical simulations in sect. 5. Section 6 is 
devoted to concluding remarks. 



2 Anisotropic lattice with improved gauge actions 

In case of improved actions that consist of plaquette and rectangular loops, 
the anisotropic lattice may be formulated in the same way as for the standard 
plaquette action |^]. The action takes the form 

59 = /?-EE^^.- + /5.EE^4., (2) 

where are plaquette and six-link rectangular loop operators in fi-v plane, 
P^, = coP(l X 1)^, + ciP(l X 2)^„ (3) 
which are constructed of the link variable Un,fj, and 

Here ^o- and Qr are the coupling constants in the spatial and temporal direc- 
tions, respectively. The action can be also written with the bare anisotropy 
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parameter as 



= E E + E E ^4.), (5) 



where f3^ = 2Njgl = and ^/^b = VdfM- 

The weak couphng hmit of the anisotropic lattice is fully discussed in 
Ref.|^. Therefore, we will summarize only equations which are necessary in 
the following studies. 

In the continuum limit — > 0{g 0), the lattice spacing and the 
coupling are related with each other by the scale parameter A through 
the renormalization group relation, 

a.A(0 = (6o5|)-'^/^''»^ exp{-l/(26o5|)}, (6) 
where bo and bi are the universal first two coefficients of /3-function, 

We calculate the effective action using the background field method [^] 
up to one-loop order, 

Seff = ^(^-^^.(0+0(5'))E4«'«- 



1 / 1 



x (8) 
+ ^[-^-CA0 + Oig')jJ2iF^, + Fl)alar. 

Effective actions with different value of the anisotropy parameter ^ corre- 
spond to different regularization scheme, but they should have the same 
continuum limit and we require A5e// = S^'jj — S\jj = 0. Then the rela- 
tions are obtained, 

1 1 



1 1 



+ (a(0-a(l)) + 0(5'), (9) 

In the following, the deviation of the one-loop quantum correction from the 
isotropic case is often employed and written as ACo-(C) = C'o-(C) — C'o-(l), 

Aa(0 = CriO - Cr{l). 
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Perturbatively, all fundamental parameters on the anisotropic lattice 
are given in terms of Ca{(,) and Cr(^). The A parameter on the anisotropic 
lattice is given by 



A(l) ' I 46„ 

The quantum correction for the anisotropy parameter r] = ^/^b is written 



as 

^2 \ S 



r/(e,/?) = ^= (^Ij' = l + yr?i(e)+0(r'), (12) 

r?i(o = a(e)-a(e), (i3) 

cwtoeiet^fiba; 2SN^y^p?^'cgiKicj^^)ciisnthffl.i:fiigiitemib^ot]^ tlfetbfi^ljG((^ 
and Co-(C) with respect to ^, 

dg-\i) dc^io dg-HO dCAO , . 

They play an important role in QCD thermodynamics ||9|, |l^, 17| as will 
be discussed in section 4.3. 



3 Perturbative calculation of Co- and Cr 
3.1 Background field method 

We calculate Co-(^) and Cr{C) in one- loop order, by applying the background 
field method. The background field method on the lattice is well known 
[18, 19 1 and therefore here we will only outline the method of the calculation 
and stress the points related to the anisotropic lattice. The gauge field is 
decomposed into a quantum field and a background one which satisfies 
the classical equation of motion, 

A gauge fixing term is introduced as, 

Sg.f. = -alarY.TT(£Dfa^{n)f. (16) 

n fi 

Here 

Z?i°^a. = ^^{uSla^in + /i)C/(°^ - a,(n)). 
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The Faddeev-Popov term resulting from the gauge fixing is 



Sp.p. = -2ala,Y.J2mD^^U{n))HD(^U{n))]. (18) 

n fj, 

The total action becomes 

Stot{ati,(t>,Bfj_) = Sg{a^,B^) + Sg,f,{a^,B^) + Sf.p.{4>, B^j^), (19) 

asH^HisrSgxipihe gauge action constructed from plaquette and six-link rect- 
It is invariant under the following gauge transformation, 

Ui^^l ^ y(n)C/i°)yt(n + ^) 

a^(n) ^ V{n)af,{n)VHn) (20) 
0(n) ^ y(n)0(n)yt(n), 

where V is an element of SU{Nc)- For the calculation of the effective action 
in one-loop order we expand the total action up to second order in and 

(2) 

Bf^, which we denote as S^J . With the help of Campbell-Hausdorff 's formula 
and a relation yexp(ia) = eyip{iVaV~^)V , we split the total action into a 
classical action and a bilinear term of and (p, 

Stot = S{B^) + {a^, B^). (22) 

The term linear in B^ is missing because of the equation of motion for the 
background field. 

Thanks to the gauge invariance of the background field, it is sufficient 
to calculate the coefficients of p^PyB^By to obtain the effective action. For 
the calculation of this term, we have applied the method explained in the 
appendix of Ref.[20|, which makes the calculation much simpler. 



It is convenient to separate the action 5"^^^^ into several parts, i.e., s'^i = 
So + Sq + ■ ■ ■ + + Sp.p.- So,Sq, - ■ ■ and Sq are symbolically expressed as 
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follows [EOI 



So 






(Aa)2 


Si 


f[B,a] 


S'l 


Aa[B,a] 


S2 


[a, a]W 


S3 


[B,a] [B,a] 


S4 


f[B,[B,a]] 


s. 


[[B,a],a]W 


Se 






aW[a,W], 



(23) 



where B and a represent B^ and a^, respectively and W and / are the 
field strength tensors of background and quantum fields, respectively. A is 
a lattice derivative when we set = 1 in Eq.(|l3). Sq and S'q are the free 
part of the action, which defines gluon propagators. Si, Si and S2 terms 
correspond to three-point diagram from which we construct the one-loop 
self-energy. S3 to Sq contribute to the tadpole self-energy. Here Sq and Si 
result from gauge fixing terms. The Faddeev-Popov term is the same as 
the previous calculations except that the anisotropy parameter ^ is included 



[18|, |ig|]. 



3.2 Effects of anisotropy parameters 

By the integration over the quantum fields, we obtain the effective action 
in one-loop order. We carry out the integration in the momentum space. 
Fourier transform of the gauge and Faddeev-Popov fields are defined as 



/TT 
a^a^{ka)exp{i{n + l/2)ka)Y[ 

d[kyay) 



arr0{n 



,(f){ka) exTp{inka) Y[ 



2tt 



(24) 



(25) 



We also define Fourier transformation of the classical field in a similar man- 
ner. 

By this Fourier transformation, the anisotropy parameters are factorized 
in the action S^g}. 



Stot,iiu{0 — -^l-tuStot,nui^) 



(26) 



7 



(2) jT^ — h 

^toifiui^) s-i'^ already given by Iwasaki and Sakai[|20| on the isotropic lattice 



and Xnu are defined as 



Xfiu — 



1 1 i t 

rrr: 
fffi 



(27) 



In this way the perturbative calculation of the anisotropic lattice becomes 
very systematic and transparent. 

For example the free part of the anisotropic improved action Sq is given 
by ^ 

So = -- / ^af,{ka)Gf,uau{-ka), (28) 




^{qukf + kf} + iqakl 

irnkf + e^l 

^{1 ~ Hj}kikj 
C{1 - q'ij}hkj 

2 sin ^A;^a 

l-ci{kl + kl) {ii^u) 
0, 



(29) 



(30) 



where the term in G44 results from the gauge fixing term and stands 
for n/i=i /-TT ^^m'^m/^^- Note that off-diagonal elements of G^p vanish for 
the Wilson action. 

Propagators D^y are defined by 

< c^^{ka)ai{k' a) >= 6ij{2TT)'^6''^\ka + k'a)D^^{ka), (31) 
and they are obtained by solving the equations, 

(32) 

The Faddeev-Popov propagator is given by 



Dp.p.iko) 



(33) 
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The explicit forms of the Si to terms in Eq. ( p3D are obtained by setting 
C2 = C3 = in formulae of Ref.|2C] and by introducing the anisotropy factors 

as shown in Eqs. (p3) and (p7|). 

— — (2) 

The effective action is obtained by integrating SfJ{a^, 4>, B^) over a^a^ 

and Faddeev-Popov fields Oq-c/), 

Then Ca{Cj and Ct(^) in Eq.(^) are obtained as coefficients of Ffj and Ff^ 
in Seff.-, respectively. 



3.3 Infrared divergence 

We will discuss here a subtle point concerning the cancellation of the infrared 
divergence. The contributions for Ca and Cr from the self-energy type 
diagram of the term < (5'i + 52)^ > have the infrared divergence. However in 
the difference given by Eqs. (^) and ([lO|) they are canceled. But numerically 
the calculation of the divergent integral is very delicate problem. In the 
numerical evaluation, we discretize the momentum integration / d^fc, and 
the divergence comes from the segment including A; = 0. We should not take 
the difference of integration with different ^ directly, because their measures 
of segment are different for isotropic and anisotropic lattices. We use the 
following method for the calculation of the Eqs. (^) and (|lO|), 

+ (Cf«"^-(l) -ci'^p-(i)). ^ ' 



where (7^'"^- and C'^^"'"''^' are the coefficients with improved and standard 
actions, respectively. 

In the first term of r.h.s of Eq.®, Ci™P (0 and C^tand.^^^ ^^^^ 
same infrared divergence and they are canceled exactly by each other. The 
second term of Eq.(|35|) can be calculated by the analytic integration of the 
Ath component of the loop momentum [^. The results do not include the 
infrared divergence and the numerical integration is stable. The divergence 
in the last term has been already calculated in Ref.|2C]. We have checked 
that our calculations for the second and the third terms of Eq. (^) coincide 
with those for Wilson and Symanzik case respectively given in Ref.g and 
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In this way the difference C^'^P'{$,) — C^'^^'{1) is calculated in numerically 
stable manner. Similar calculation has been done for Cl^^'{^) — C^™'^'(1). 

4 Results of one-loop calculation 



Values of Ca{Cj — Ca{'^) and Cr(C) — Cr(l) are given in Table |l|for Symanzik, 
Iwasaki and DBW2. ^ is varied from 1 to 6, since these anisotropy parame- 
ters are often used in Monte Carlo simulations on anisotropic lattices. 







Symanzik action 


Iwasaki action 


DBW2 action 






AC 


Aa 


AC 


Aa 


AC 


Aa 




1.0 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 




1.1 


0.00541 


-0.00092 


-0.00201 


0.00614 


-0.01124 


0.01464 




1.5 


0.02072 


-0.00173 


-0.00908 


0.02327 


-0.04875 


0.05150 


SU(2) 


2.0 


0.03180 


-0.00085 


-0.01564 


0.03612 


-0.08400 


0.07510 




3.0 


0.04304 


0.00103 


-0.02406 


0.05016 


-0.13384 


0.09537 




4.0 


0.04862 


0.00159 


-0.02934 


0.05613 


-0.16932 


0.09766 




5.0 


0.05168 


0.00059 


-0.03340 


0.05692 


-0.19695 


0.09055 




6.0 


0.05314 


-0.00162 


-0.03725 


0.05444 


-0.21970 


0.08000 




1.0 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 




1.1 


0.00930 


-0.00253 


-0.00427 


0.01045 


-0.02087 


0.02584 




1.5 


0.03602 


-0.00701 


-0.01852 


0.03967 


-0.08983 


0.09206 


SU(3) 


2.0 


0.05571 


-0.00829 


-0.03119 


0.06184 


-0.15389 


0.13609 




3.0 


0.07596 


-0.00856 


-0.04693 


0.08679 


-0.24348 


0.17729 




4.0 


0.08614 


-0.00951 


-0.05644 


0.09843 


-0.30658 


0.18776 




5.0 


0.09182 


-0.01216 


-0.06350 


0.10163 


-0.35533 


0.18221 




6.0 


0.09462 


-0.01625 


-0.07012 


0.09953 


-0.39532 


0.17034 



Table 1: a and a for Symanzik, Iwasaki and DBW2 actions. Here Aa = 

a(0 - Ca(i) and Aa = a(0 - a(i). 



4.1 The A ratio 

When we calculate physical quantities, we must take into account the vari- 
ation of the scale due to A(^). In weak coupling regions it is given by 

^Although the values of C'^*""'' (C) - C^*""'*!!) were not given in Ref.§, we have 
calculated them from Table and formulae in Ref. H 
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Eg. (pA|) . The A ratio is calculated as a product of three factors, 

Ajmp. iO _ Ajmp.(C) ^ ^Stand. (Q ^ ^Stand.i^) ^gg-j 
A/mp.(l) ^StandXO ^StandX^) ^Imp.{^) 

The numerical results are given in the Table |2| for Symanzik, Iwasaki and 
DBW2 and are shown in Fig. 1. 







A/„p.(0/A/mp.(l) 




A/mp.(0/A/mp.(l) 






Symanzik 


Iwasaki 


DBW2 




Symanzik 


Iwasaki 


DBW2 




1.0 


1.00000 


1.00000 


1.00000 




1.00000 


1.00000 


1.00000 




1.1 


0.97614 


0.97803 


0.98187 




0.97599 


0.97809 


0.98234 




1.5 


0.90282 


0.92641 


0.98529 




0.90111 


0.92691 


0.99204 


SU(2) 


2.0 


0.84653 


0.89562 


1.04907 


SU(3) 


0.84350 


0.89583 


1.06597 




3.0 


0.78873 


0.86890 


1.23013 




0.78512 


0.86670 


1.26812 




4.0 


0.76309 


0.86566 


1.47071 




0.75955 


0.86010 


1.53177 




5.0 


0.75471 


0.88107 


1.77320 




0.75133 


0.87209 


1.86138 




6.0 


0.75778 


0.91159 


2.12140 




0.75483 


0.89981 


2.24211 



Table 2: A parameter ratio for improved actions. 



2 



0.5 



• — • standard Action 








■ ■ Symanzik Action 








♦ - - ♦ lwasal<i Action 








A- - A DBW2 Action 




y 

y 















































♦ 


— •♦ 


*> 




■ 


• 


- !! 






■ 





1 2 3 4 5 6 

4 



Figure 1: The ratio of A parameter for SU(3). 
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4.2 The 7] parameter 

The rj parameter is given by Eqs. ([l2|) and (|l3|). Because of the gauge 
invariance we can see Ca{l) = Ct-(1) at the isotropic case, and therefore these 
coefficients do not appear in the definition of rj. But to cancel the infrared 
divergence we calculate Eq.(13) asr?i(C) = {C^{£,)-C„{l))-{Cr{0-Cr{l)). 
We plot r]i as a function of ci in Fig. 







Symanzik 


Iwasaki 


DBW2 




Symanzik 


Iwasaki 


DBW2 




1.0 


0.00000 


0.00000 


0.00000 




0.00000 


0.00000 


0.00000 




1.1 


0.00633 


-0.00816 


-0.02588 




0.01183 


-0.01472 


-0.04671 




1.5 


0.02246 


-0.03235 


-0.10025 




0.04303 


-0.05820 


-0.18190 


SU(2) 


2.0 


0.03266 


-0.05177 


-0.15911 


SU(3) 


0.06400 


-0.09304 


-0.28998 




3.0 


0.04201 


-0.07423 


-0.22921 




0.08452 


-0.13372 


-0.42077 




4.0 


0.04702 


-0.08548 


-0.26698 




0.09565 


-0.15487 


-0.49434 




5.0 


0.05109 


-0.09032 


-0.28751 




0.10398 


-0.16513 


-0.53755 




6.0 


0.05476 


-0.09170 


-0.29971 




0.11088 


-0.16966 


-0.56567 



Table 3: The r]i for improved actions. 



The parameter rji changes sign around ci = —0.18 ~ —0.19. This means 
that in weak coupling regions, there is no renormalization for the anisotropy 
parameter ^ for this action. The interesting point is that the (3 dependence 
of the T] with Iwasaki and DBW2 actions is opposite to those with standard 
and Symanzik actions; As (5 decreases t] decreases for Iwasaki and DBW2 
action, while it increases for standard and Symanzik actions. This is a new 
feature and we shall confirm it non-perturbatively in the next section. 



4.3 Anisotropy coefRcients 

The anisotropy coefficients, which are the derivatives of spatial and temporal 
gauge couplings with respect to the anisotropy ^, are calculated as, 

= -^{Ci^^P-iO - '^"'^■(0) + ^(Cf - Cf '^"'^■(1)). (37) 

In this manner, we are free from the infrared divergence, and the numerical 
evaluation is stable as in Sect. 3.3. 
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Figure 2: The r]i as a function of ci for SU(3). 



Prom the invariance of the string tension on the isotropic and anisotropic 
lattice, Karsch has derived the following sum ruleP], 



UN, 
487r2 



(C = l,a^O). 



(38) 



The same arguments are applied to improved actions. We show results in 
Table § and Fig. ^. This sum rule is satisfied quite well. 

An interesting point is that individual terms dg~'^ /d^ and dg~'^ /d^ also 
change sign as — ci increases. For Iwasaki and DBW2 actions, dg~'^/d^ 
becomes positive while it is negative for standard and Symanzik actions. 
These anisotropy coefficients have a contribution to QCD thermodynam- 
ics. For example, Okamoto et al. [p2|| used our perturbation results to study 
the energy and pressure with Iwasaki action and had no negative pressure 



problem, contrary to Wilson action case|23] 
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5 Numerical results in weak coupling regions 



In the previous section we have found by the perturbative calculation that 
the ratio of the renormalized and bare anisotropy, r], becomes less than one 
for Iwasaki and DBW2 actions. 

Since the r/ is important in QCD simulations on the anisotropic lattice, 
we will study its behavior further by numerical simulations. 

Numerically the r] parameter is calculated from the relation |2^ , 

V = C/Cb. (39) 

The anisotropy S^b appears in the action given by Eq.(^) while the renor- 
malized anisotropy ^ is defined by 

e = ajar. (40) 

For the probe of the scale in the space and temperature direction, we use 
the lattice potential in these directions, which is defined by 

where Wari^yt) is the Wilson loop of the size / x t in the temporal plane. 
Similar formula holds for the potential in space direction. We fix ^ = 2, and 
calculate the ratio at a few S,b points, 



Then we search for the point i? = 1 by interpolating ^_b[0, 24|. In Ref.|25], 
an extensive study is done for the determination of S,b- 





SU(2) 


SU(3) 








d9jmk=i 




Symanzik 
Iwasaki 
DBW2 


0.058406 
-0.020569 
-0.117103 


-0.01196 
0.067009 
0.163543 


0.10025 
-0.044305 
-0.21793 


-0.030365 
0.113965 
0.287592 



Table 4: The anisotropy coefficients for improved actions. The results are 
checked by estimating the derivative numerically by using the results with 
^ = 1.05,1.1. 



14 



-0.2 



dg;'/d^ 
xdg;'/d^ + dg;'/d^ 



-0.4 



-1.5 



-0.5 



Figure 3: Anisotropy coefficients with the standard and improved actions 
for SU(3). Dotted hne represents the r.h.s of Eg . (|3^) , which is about 0.0697 
for SU(3) case. The cross symbol stands for the l.h.s of Eq.(^). 



The simulations are done on the 12^ x 24 lattice. Numerical results with 
large /? region( /3 > 10) together with perturbative ones are shown in Fig. 
^. They agree with each other at large f3 region. The tj parameter decreases 
as P decreases for Iwasaki and DBW2 while it increases for standard and 
Symanzik actions. 



6 Concluding remarks 

We have calculated the QCD scale parameter A, which is shown in Table § 
and Fig. ||. Ratios A(^)/A(l) for standard, Symanzik and Iwasaki actions 
have very similar behavior; they are slightly less than one, but the behavior 
of A ratio of the DBW2 is quite different from those of other actions. DBW2 



action is expected to be very near to the renormalized trajectory 1 15], and 
may have a special feature. It may be interesting to study the ci dependence 
of the A parameters between Iwasaki and DBW2 actions in more detail. 
The r] parameters and anisotropy coefficients, the derivatives of the cou- 
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Figure 4: Perturbative and non-perturbative results of r/ as a function of 
p. Data for Wilson actions are taken from Ref.[24| 



pling constants with respect to the anisotropy parameter ^ are calculated 
in one-loop order for improved actions. We have found that the rji in the 
Eq.([T3|) changes sign for ci ~ —0.18 ~ —0.19 as shown in Fig. ^. The 
is positive for the standard plaquette and Symanzik actions while it be- 
comes negative for Iwasaki and DBW2 actions. This behavior is confirmed 
non-perturbatively by the numerical simulations in weak coupling regions 
as shown in Fig. |^. 

We have also found that the anisotropy coefficients change sign as — ci 
increases. For Iwasaki and DBW2 actions, dg~'^/d(, is positive while for 
standard and Symanzik actions it is negative. These may be good properties 
for the thermodynamics using those improved actions |p2|. 

We have found that rj obtained by the perturbation calculations is close 
to one in the region ci ~ —0.18 ~ —0.19 for ^ = 1 ~ 6 0. A natural question 
is whether this is true also at intermediate and strong coupling regions. 



^The value is sligtly shifted from that in Ref. After lattice'99, we have done an 
extensive study of the numerical calculation of the loop integral, so that we can safely 
extrapolate to the continuum integral. 
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Parts of the results on the rj parameter are reported at lattice '99 at Pisa|21]. 
The study of the lattice spacing a on the anisotropic lattice in intermediate 
coupling regions has been started. Detailed results will be reported in the 
forthcoming paper. Moreover, with these fundamental properties of the 
improved actions on anisotropic lattices, we are going to start simulations 
of heavy quark spectroscopy, transport coefficients of quark gluon plasma 
etc. on these lattices. 
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Appendix 

In this appendix we show all the data used in Eq. (|35| ) . In Table |5| the data 
of the first term in Eq.([35|) are summarized, and in Table |^ we report the 
data for the standard action. The values of the third terms in Eq.(^5|) are 
given by isotropy cases in Table ||. 
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Symanzik action 


Iwasaki action 


DBW2 action 




? 




5Cr iO 








6Cr{0 




i.U 


-0.13173 


-0.13173 


i~\ o r\ <^ n 

-0.32668 


-0.32668 


-0.57505 


-0.57505 




i.i 


-0.13690 


-0.12672 


-0.33928 


-0.31461 


-0.59687 


-0.55447 




i.o 


-0.15149 


-0.11355 


-0.37625 


-0.28349 


-0.66429 


-0.50362 


CTT/0\ 

bU (2) 


2.0 


-0.16130 


-0.10393 


-0.40371 


-0.26190 


-0.72042 


-0.47128 




o.U 


-u.ioyuo 


-U.UyzDi 


-U.4oii0 


n OQQ/1 /I 
-U.Zoo44 


-U. loyZy 


n /I /1 1 f^r\ 






-0.17135 


-0.08643 


-0.44427 


-0.22685 


-0.83262 


-0.43368 




5.0 


-0.17218 


-0.08359 


-0.45223 


-0.22221 


-0.86414 


-0.43694 




6.0 


-0.17292 


-0.08300 


-0.45828 


-0.22188 


-0.88909 


-0.44468 




1.0 


-0.23211 


-0.23211 


-0.57349 


-0.57349 


-0.99987 


-0.99987 




1.1 


-0.24156 


-0.22293 


-0.59652 


-0.55132 


-1.03950 


-0.96232 




1.5 


-0.26836 


-0.19832 


-0.66429 


-0.49302 


-1.16198 


-0.86701 


SU(3) 


2.0 


-0.28649 


-0.17988 


-0.71477 


-0.45112 


-1.26386 


-0.80326 




3.0 


-0.30099 


-0.15771 


-0.76527 


-0.40373 


-1.38820 


-0.73962 




4.0 


-0.30524 


-0.14509 


-0.78921 


-0.37852 


-1.46573 


-0.71558 




5.0 


-0.30674 


-0.13848 


-0.80344 


-0.36606 


-1.52166 


-0.71187 




6.0 


-0.30799 


-0.13586 


-0.81412 


-0.36144 


-1.56570 


-0.71702 



Table 5: and Cr for Symanzik, Iwasaki and DBW2 actions. Here 
5C^{^) = Cl^'P-iO - C|*"'^'^-(e) and 6Cr{0 = - ""'^■(C). 
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